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The program was conducted by TRW Systems Group of TRW, Inc., 

Redondo Beach, California. Dr. Bruce D. Marcus served as Program 
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1.0 INTRODUCTION 

For the last several years TRW Systems Group, under contract 
to NASA-ARC, has performed an extensive research and development program 
in variable-conductance heat-pipe technology. The treatment has been 
comprehensive, involving theoretical and/or experimental studies in 
hydrostatics, hydrodynamics, heat transfer into and out of the pipe, fluid 
selection, and materials compatibility, in addition to the principal 
subject of variable-conductance control techniques. Efforts were not 
limited to analytical work and laboratory experimentation, but extended 
by the development, fabrication and test of spacecraft hardware, highlighted 
in the successful flight of the Ames Heat Pipe Experiment on the OAO-C 
spacecraft. 

Most of the program's accomplishments have been previously 
documented in a series of reports and publications. Early theoretical and 
design developments appear in References [1-1, 1-2, 1-3, 1-4, 1-5]. 

Later fundamental work was published in References [1-6, 1-7, 1-8]. Hardware 
development and application efforts were documented in References [1-9, 

1-10, 1-113, and a computer program for desi ciing and predicting performance 
of gas loaded heat pipes was presented in Reference [1-12]. 

This document represents the final report on the contract 
and presents the results of TRW's latest efforts. It does not attempt 
to present all of the work accomplished on the program, since most has 
been previously published. However, to provide a summary of the program's 
scope and a guide as to the location of previously published information, 
copies of the Tables of Contents of all earlier reports are included in 
the Appendix. 

The particular studies that are covered in this report fall 
into four areas as follows: 

1) An experimental and theoretical study was made 
of gas generation in methanol/stainless-steel 
heat pipes for the purpose of establishing a 
scaling law for the accelerated testing of such 
heat pipes (Section 2.0). 
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2) The TRto GASPIPE computer program, previously ? 

developed on this contract (Ref. [1-12], was 

extended to include two condenser sections 
with different properties, which enables the 
program to be easily applied to gas-loaded 
heat pipes having a secondary short condenser 
(cold trap) adjacent to the gas reservoir as 
well as an adiabatic section (Section 3). 

3) Theoretical and experimental work was carried \ 

$ 

out on the effect of noncondensable gas on 

arterial performance. Pressure fluctuations in \ 

gas-loaded ammonia heat pipes were studied, and I 

a glass heat pipe was fabricated to study a new j 

method of venting noncondensable gas during prim- f 

ing (Section 4.0) . 

4) Two research heat pipes (and two spares) were 
designed, fabricated and tested for a forth- 
coming 6SFC sounding-rocket experiment to study 
arterial priming in zero gravity. These heat 
pipes have internal thermistor instrumentation 
that provides information on the priming process 
(Section 5.0). 

: 
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2.0 DEVELOPMENT OF GAS GENERATION SCALING LAWS FOR METHANOL/ 

STAINLESS-STEEL HEAT PIPES 

Heat pipes are rapidly becoming a serious design element 
ir. the solution of many spacecraft thermal control problems where long 
periods of trouble free performance are reouired. For design purposes, 
therefore, there is interest in determining methods for estimating the 
operating lifetime. What is meant by "operating lifetime" depends on the 
type of heat pipe and the requirements placed on it in a specific appli- 
cation. With high temperature liquid metal heat pipes, for example, 
structural failure sometimes occur and this is a well-defined termination 
of life. In low temperature heat pipes applicable to spacecraft thermal 
control, catastrophic failures rarely occur. Instead, heat-^ice performance 
continuously degrades as a result of (1) chemical reaction or decomposition 
of the working fluid with the generation of noncondensable gas, or (2) 
corrosion and erosion of the container and wick. 

In an ordinary heat pipe all noncondensable gas is swept to 
the condenser end, forming a diffusion barrier to vapor flow and effectively 
reducing the available condenser area. In gas controlled, variable conduc- 
tance heat pipes, the generation of additional noncondensable gas raises 
the operating temperature of the heat pipe above design conditions. Similar 
effects can result from a change in the chemical composition of the working 
fluid by virtue of a change in its vapor pressure as a function of temperature. 

Corrosion and erosion of the container and wick can be manifested 
as a change in the wetting angle of the working fluid as well as the 
permeability, porosity, or capillary pore size of the wick. Solid 
precipitates resulting from corrosion and erosion are transported by the 
flowing liquid to the evaporator region where they are deposited when 
the liquid vaporizes. This leads to increased resistance to fluid flow 
in the evaporator, resulting in a decrease in the heat transport capacity 
of the heat pipe. 

With these failure mechanisms, where continual degradation 
occurs, the operating lifetime ca 1 ” be defined as that period of time beyond 
which the operation of the heat pipe is below design specifications. Some 
heat pipe laboratories have been performing "life tests" under which heat 
pipes are held at normal operating conditions for many thousands of hours 
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to determine the "operating lifetimes". This approach has limited 
applicability, however, for heat pipes which are required to function 
well for long periods of time. Progress in heat-pipe development will 
be impeded if each time a new material combination, fabrication technique, 
or cleaning procedure is used, life tests of 10,000 hours or more are 
required. There exists, therefore, a major impetus for understanding 
the chemical and corrosion mechanisms, and developing scaling laws 
for the life-limiting processes. The need for achieving very long lived 
high reliability heat-pipe systems in a rapidly charging technology 
necessitates the employment of accelerated testing techniques. 

This program was the second in a continuing effort to under- 
stand compatibility problems as they relate to heat-pipe operating 
lifetimes. The initial study [2-1, 2-2] in this series was carried out 
using nickel/water heat pipes as a feasibility demonstration to assess 
the difficulties in establishing a scaling law relationship for a 
specific example. From a study of hydrogen evolution in nickel/water 
heat pipes under accelerated (high) and reference (low) operating 
conditions it was found possible to accurately predict, with a scaling 
law, the lifetime cf a heat pipe operating at reference conditions from 
data taken at accelerated conditions. It was also found that the same 
form of scaling law correlated the data of Petrick [2-3] on stainless- 
steel/water heat pipes. These results suggested that the formulation of 
crying laws for accelerated testing of life-limiting processes is indeed 
feasible and that similar methoGS might be applicable to other heat-pipe 
systems as wel 1 . 

The present program is concerned with developing scaling laws 
for gas generation in stainless-steel/methanol heat pipes. Previous life 
testing of such heat pipes has shown this combination to be compatible 
up to 145°F (63°C) , but gas generation has been observed in these heat 
pipes at higher temperatures. As stainless-steel /methanol heat pipes 
represent the most desirable combination for many variable conductance 
heat pipes, it was important that a study be made of the gas generation 
behavior at higher temperatures. Before this program was carried out, 
even the high temperature limit of operation, before the initiation of 
significant gas generation, was unknown. 
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The primary objective was to formulate an accurate scaling law, 
based on a chemical and corrosion model of gas evolution mechanisms, for 
stainless-steel /methanol heat pipes, which predicts the usable lifetime 
at reference (low) operating conditions from data taken at accelerated 
(high) operating conditions. This method of accelerated life testing 
is based on extrapolation from accelerated conditions to predict behavior 
at reference (normal) conditions. It is thus explicitly ass >t the 

chemical and physical mechanisms responsible for heat-pipe ■ gradat 
at accelerated conditions are the same in nature at the reference conditions. 
Similar ideas are used for other types of accelerated corrosion testing 
[2-4]. It is emphasized that this work was not concerned with achieving 
compatibility between stainless steel and methanol. However, strict 
control was maintained over the starting materials and fabrication 
techniques, not to eliminate gas generation, but to allow the separation 
of temperature and fluid circulation effects. 

2.1 Accelerated Life Testing 

2.1.1 Heat Pi k e Materials and Fabrication 

Sixteen heat pipes were fabricated using 304 stainless-steel 

materials. Fabrication procedures which are standard at TRW for construction 

of stain! ess-steel /methanol heat pipes were used in orcer that the results 

of accelerated life testing would be directly applicable to actual heat 

pipes built by TRW. Container tubes were 19.5" in length with 1/2" O.D. 

and 0.020" walls. Two layers of 150 mesh (2.6 mil diameter) screen were 

installed and the end cops were machined from 1/2" O.D. rods. Closure 

tubes were cut from 1/8" O.D., 0.020" wall tubing. Care was taker, to 

ensure that all the heat pipes were as nearly t. e same as possible in 

terms of materials and construction procedures. All containers were cut 

from 304 stainless steel tube of the same heat number and the same was 

true of the closure tubing and rods from which the end caps were machined. 

Screens were cut from a single sheet of wire mesh. Welding was performed 

with argon in the tube using 308 stainless-steel welding rod and care 

was taken to use the same temperature and coi te the weld in the same 

length of time for each heat pipe. After the screens were installed and 

the parts cleaneG, one eno cap was welded, and in this condition all the 

-4 

parts were vacuum fired for 1 hour at 100C°C below 10 torr. The remaining 
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welds were then completed and the heat pipes were filled with 5.5 ml of 
methanol (Matheson Coleman & Bell, spectroqual ity grade containing a 
maximum of O.C5% water) and instrumented with 6" heaters. A schematic 
diagram of a heat pipe prepared for accelerated testing is shown in 
Figure 2-1. After filling, one heat pipe was found to be unusable due 
to excess gas in the pipe, probably resulting from a leak during the 
filling procedure. All other heat pipes were found to contain initially 
only a small amount of gas. Eleven copper-constantan thermocouples were 
placed at 3/4" intervals along the condenser, and one thermocouple was 
placed in the adiabatic section, for the purpose of measuring the non- 
condensable gas content by means of the temperature profile. An additional 
thermocouple was placed on the end cap at the heater end to test for burn- 
out. Insulation consisted of an inner layer of J-M micro-fibers felt and 
an outer layer of J-M aerotube. 

2.1.2 Measurement of Noncondensable Gas Evolution 

As noncondensable gas is evolved during operation of a heat 
pipe, it is carried to the condenser end causing blockage and a consequent 
temperature profile along the wall. The amount of gas present may be 
calculated from the temperature profile assuming ideal gas behavior. If 
the condenser end is divided into N equal intervals and the temperature 
at the center of each interval is , then under steady-state cond?*‘~ . 
the number of lb moles of gas n is given by the ideal gas law as: 


N 


aV 

n = U- z 

R 1 = 1 '1 


^ ■ 


(2-1) 


where aV is the volume of each interval, R is the gas constant, and 


P gi ’ P ,a * P ,1 


(Z-Z) 


th 


is the partial pressure of gas at the center of the i ' interval. In 
(2-2), P ya Is the total pressure (the vapor pressure corresponding to 
the temperature in the adiabatic section) and \i is the vapor pressure 
in the i th interval. 


A computer program was used to calculate the quantity of gas 
In a heat pipe at any given time from the measured steady-state wall 
temperature profile. This method Is based on the assumption that the 

6 




FIGURE 21. Schematic Diagram of Heat Pipes Subjected co Accelerated Life Testing. 
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wick surface temperature, and hence the vapor-gas mixture, is very close 
to the wall temperature in the gas-b1ock*a region of the condenser. This 
has been found to be a valid assumptio . from the study of gas controlled 
pipes [2-5]. 

In practice, each pipe was divided into 3/4" elements with a 
thermocouple placed in the center of each interval. The thermocouple 
temperature readings were input directly into the computer program, 
which carried out the operations indicated by Eqs. (2-1) and (2-2) and 
printed out the total number of lb moles of gas in the pipe. Trial 
calculations have indicated that the discrepancy between using a 0.5" 
element and a 1.0" element was less than one percent. 

2.1.3 Results of Accelerated Testing 

Based on previous studies of gas evolution in heat pipes 
[2-1 rhrough 2-3, 2-6, 2-7], it was assumed that the gas generation rate 
would be a strong function of tne operating (vapor) temperature. For 
this reason all heat pipes were tested in a constant temperature chamber. 
The temperature of the chamber was unaffected by convective air currents 
in the room and could be held at -0.5°C of the set point. 

Initial accelerated life testing was begun with three heat pipes 
operated at 66°C, 9J°C, and I21°C. Other than a small amount of gas which 
was present from the filling procedure, no detectable additional amount 
was generated after 39 days of operation. Initial testing was then 
continued with the temperature increased to 177°C, 204°C and 232°C. 
Temperature profiles were readily generated at these temperatures over 
a two week period, the gas evolution rate increasing with increasing 
temperature. 

During initial testing, it was found that burnout occurred at 
temperatures above 204°C even with the heat pipe in a vertical reflux 
position. Thus data taken with the heat pipe operated at 232°C was not 
usable. Calculations of the vapor mass as a function of temperature 
indicate that burnout above 204°C resulted from liquid depletion. With 
the heat pipes operating at 177°C and 204°C, the tilt (in the reflux 
position) was adjusted to the point just above the angle which produced 
the start of burnout, as evidenced by a rise in the end-cap thermocouple 
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temperature. This procedure was used with all subsequent heat pipes to 
ensure maximum fluid flow in the wicks. 

After completing measurements on the initial set of heat pipes, 

which determined the temperature range of interest, a set of four heat 

pipes were operated (in a reflux position) at essentially the same 

temperature (179.2 - 179. 3°C) but at different flow rates (11.0, 12.4, 

16.4 and 19.6 watts), to investigate the possibility of a flow rate 

dependence in the gas generation rate. For a given power, the condenser 

length was adjusted to achieve the desired temperature. The results are 

shown in Figure 2-2, where the curves are drawn to a least-squares fit 

to the data (as with all subsequent graphs). The gas generation data 

indicates an initial parabolic time dependence, (time)^, followed by a 

-9 

linear time dependence after the generation of approximately 6 x 10 lb 
moles of gas. Although the power level was varied by nearly a factor of 
2, no flow rate dependence was apparent in the data. Thus, any high 
temperature flow rate dependence must be small, particularly when compared 
to the temperature dependence. 

Having completed the investigation of the flow rate dependence, 
the program was continued with a study of the time and temperature depen- 
dence. Because the initial studies had indicated that the gas evolution 
rate is very low at low temperatures it was apparent that the investi- 
gation of heat pipes at truly reference (normal) operating conditions 
did not fit within the time frame of this program. Thus all the 8 
remaining heat pipes were studied at a series of high (accelerated) 
temperatures. One of the 179.3°C heat pipes studied for flow rate 
dependence was included, giving 9 heat pipes at the following series 
of temperatures: 120.0, 130.0, 140.0, 149.9, 160.0, 170.2, 179.3, 

189.7, and 198. 2°C. The indicated temperatures are average values of 
recorded temperatures taken over the entire exposure period at approximately 
24 hour intervals. No temperature variations greater than - 0.5°C were 
bserved ove. any significant period of time. Only the three lowest 
temperature heat pipes could be operated in a heat pipe mode without 
burnoi't. The remainder were operated In a reflux position. A 3/4" 
condenser was exposed in each case giving corresponding flow rates of 
7.0, 7.7, 8.9, 10.0, 10.7, 11.4, 12.6, 13.6 and 15.9 watts. 
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FIGURE 2-2. Gas evolution in methanol/stainless-steel heat pipes operated 
at different power levels to investigate the possibility of a 
flow rate dependence. 
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The gas generation curves are shown in Figures 2-3 and 2-4. 

As with the heat pipes studied for flow rate dependence, the gas generation 
data appears to obey a parabolic time dependence up to approximately 

_q 

6 x 10 lb mole of gas and a linear time dependence beyond this point. 
Thus, the data points in Figure 2-4 are fitted to a parabolic curve, 
while the higher temperature data in Figure 2-3 are fitted to a linear 

_g 

time dependence above 6 x 10 lb mole and a parabolic dependence below. 

The data from the 120°C heat pipe is not included because its behavior 
was significantly different than the other heat pipes. It showed an 
abnormally high gas generation rate, possibly resulting from a very small 
leak or from impurities not contained in the other heat pipes. 

The scatter in uie data appears to be greater for the lower 

temperature pipes in which only a small amount of gas was generated. 

This may result from the fact that the results of the computer program 

are quite sensitive to small variations in the thermocouple readings, 

-9 

particularly for small amounts of gas (10 lb mole range). Thus the 
parabolic least squares fits in the range of low gas content were not 
as good (lower percent of determination) as the linear fits above 
6 x 10‘ 9 lb mole. 

2.2 Gas Generation Model and Analysis 

Stainless steel is reported to undergo uniform corrosion in 
methanol at a sma 1 ! linear rate (less than 0.1 mpy below 118°C) aftw. 
long duration exposure [2-8]. No studies of the initial corrosion 
behavior or of the temperature dependence of the corrosion rate could 
be found in the literature. In the presence of a catalyst, methanol 
undergoes dehydrogenation to form formaldehyde gas and hydrogen gas 
[2-9]: 

CH 3 0H -+ HCHO + H 2 . (2-3) 

If oxygen is also present, the overall reaction yields formaldehyde 
and water [2-10]: 

CH 3 0H + 1/2 0 2 - HCHO + H 2 0, (2-4) 

and oxidation of formaldehyde produces formic acid: 

HCHO + 1/2 0 2 -*■ HCOOH. (2-5) 
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FIGURE 2-4. Gas generation in methanoV^tainl ess-steel heat pipes showing temperature 
dependence of gas evolution in the parabolic region. 
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Methanol/stainless-steel heat pipes might be expected to generate 
formaldehyde and hydrogen as noncondensahle gases, and perhaps formic 
acid as the reaction proceeds. In fact, an analysis of the gas generated 
in methanol /stain! ess-steel heat pipes life tested for 5000 hours between 
42-55°C by kreeb [2-11] showed the presence of formaldehyde plus minor 
amounts of H 2 , CH^, and CO. 

The corrosion rate of stainless steel in formaldehyde is 
slightly greater (less than 2 mpy at 24°C) [2-12] than in methanol, 
but the rate increases by more than an order of magnitude in the presence 
of formic acid (20-50 mpy at 52 6 C) [2-12]. Uniform corrosion and pitting 
occur with both of these corrosion reactions [2-8]. Micrographs of 
sections of methanol /stainless- steel heat pipes life tested at 57°C 
for 5500 hours by Groll et al. [2-13] showed the presence of a surface 
film (uniform corrosion) and grain boundary attack. A phenomenological 
gas generation model is considered below which incorporates corrosion and 
oxidation theory and contains parameters which may be determined by 
experiment. 

It is assumed that uniform corrosion occurs at all stainless- 
steel surfaces by the operation of a great number of microscopic galvanic 
cells [2-14]. The results of this study indicate that passivating film 
growth with a parabolic time dependence probably occurs initially from 
the corrosion of the stainless-steel surface in the presence of methanol 
accompanied by the evolution of a proportional amount of noncondensable gas. 
Once the passivating film is formed, this corrosion product may then act 
as a catalyst for reactions (2-3) and (2-4), giving a linear time dependence. 
Reaction (2-4) may predominate, with the oxygen coming from the small water 
content of the methanol , based on the results of Keeb [2-11] who found 
formaldehyde to be the main constituent in the gas. Reaction (2-5) 
probably plays only a minor role in the corrosion processes. 

The transition from the passivating region to the catalytic region 
should depend on the temperature and the amount of corrosion product present 
per unit area (effective film thickness) or, equivalently, on the quantity 
n c /A, where n c is the critical gas content and A is the Internal corroded 
surface area. Since the film of corrosion products may be porous and 
discontinuous, the appropriate area A may be difficult to define. Thus, 
a critical time -t c to the transition point can also be defined, as discussed 
later, which may be a preferable criterion for the transition. 


14 



1 311 1 -6060-RU-00 


Generally accepted theories of film growth d-ring oxidation of 
metal surfaces assume diffusion through the film of cations away from the 
metal surface and anions toward the metal surface [2-15, 2-16]. The 
diffusing ions migrate from one position of minimum potential energy to 
the next. If Q is the height of the barrier between two potential energy 
minimums, the probability that an Ion will pass over the barrier is 
proportional to e -( ^ kT , where Q is called the activation energy, k is 
Boltzmann's constant, and T Is the absolute temperature. This is the 
temperature dependence found experimentally In solid state diffusion 
[2-17]. It should be mentioned that the simple Interpretation of Q 
given above is only one of several physical interpretations, which depend 
on the particulars of film growth theory. 


Oxidation theory [2-15] predicts passivating film growth will 
occur with a parabolic time dependence and an exponential temperature 
dependence. Assuming a proportional amount of gas is evolved in the 
process, gas generation in the passivating region. t <t c , should 

be given by 


n(t,T) = B 1 At 1/2 e‘ Q l /kT , (2-6) 

where n is the number of lb moles of gas, t is the time, A is the total 

p 

internal area of stainless steel in contact with the methanol (166.6 in for these 
heat pioes),and B-j is a constant characteristic of the corrosion process. This time 
dependence has been found to describe hydrogen evolution from steel in boiling 
water [2-18] but is not common to the (long duration) corrosion of metals 
generally. Various other forms of time dependencies are also predicted 
by theory depending on the particular assumptions made [2-19, 2-20]. Over 
long exposure periods, uniform corrosion with a linear time dependence 
occurs more commonly. 

According to the gas generation model under consideration, after 
sufficient corrosion product has been produced to act as a catalyst for 
reactions (2-3) and (2-4), equivalent to a quantity of gas n c /A lb mole per unit 
area at the transition point, formaldehyde and hydrogen are produced. These 
catalytic reactions are considered to predominate beyond the transition point, 

and would be expected to obey a linear time dependence with a characteristic 
activation energy Qg* giving 
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- % 


n(t,T) * B 2 Ate~ Q 2 /kT , (2-7) 

for n/A>-^-, t >t c , where Is a constant characteristic of the catalytic 
reaction. Here Qg represents the potential barrier between reacting molecules 
on the catalytic surface. This temperature dependence is characteristic of 
many physical and chemical reactions [2-15]. The area A may actually be 
somewhat larqer than geometrical Internal area depending on the surface 
roughness of the stainless steel or corrosion product surfaces . 


With dry corrosion the activation energies are generally larger 
than with wet corrosion, as Is apparent from the activation energies 
shown in lable 2-1. An explanation for this may be that the effective 
potential barrier is lowered by the electric field across the film 
created by the local corrosion cell. 


The use of this gas generation model in accelerated life 

testing is that the parameters B and Q con he determined experimentally 

from data taken under accelerated conditions by plotting log 9n vs. 

dt 

l/T. Having determined these parameters by measuring the gas 
evolution at accelerated conditions the gas evolution at any time can 
be calculated from (2-6) and (2-7) ror heat pipes operated under normal 
conditions. The quantities t c and n £ /A can be calculated from (2-8) and 
(2-9), as discussed below. 


Least squares fits to the parabolic data are plotted n vs t 1/2 

in Figures 2-5 through 2-7. Plotting log II > vs ] /T results in the 

3t'^ z 

curve shown in Figure 2-8, indicating gas generation In the passivating 
region is described by (2-6), within the accuracy of the data. Calculating 
the parameters and B 1 from the slope and Intercept, respectively, 
results in 


Q 1 = 6.03 x 10’ 20 joules, 

B-j s 3.66 x 10”® lb moles/hr^ 2 in 2 . 
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Table 2-1. Activation energies for corrosion in gaseous and liquid environments. 
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♦Calculated from data contained in referenced papers. 
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FIGURE 2-5. Gas generation In me thanol/stalnl ess-steel heat pipes 
showing temperature dependence In parabolic region. 
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FIGURE 2-6. Gas generation In methano]/sta1nless-steel heat pipes 
In parabolic region. 
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3n 1 

A plot of log ^ vs /T In the linear region results in the 
curve shown in Figure 2-9, showing good agreement with (2-7). Calculating 
the parameters for the gas generation behavior beyond the passivating 
region results in 


Q 2 = 26.7 x 10' 20 joules, 

B 2 s 3.79 x 10 5 lb mole/hr in 2 . 


The transition point t c between the passivating and catalytic 
regions may be defined as the point beyond which catalytic reactions pre- 
dominate, i.e., the time at which the gas generation rates, |r-, are equal. 
Equating the slopes, of Eqs. (2-6) and (2-7) results in: 



exp 


‘2(Q 2 - Q.,) 

. — ET 


(2-8) 


for the critical time. The critical gas content per unit area (equivalent 
to a critical film thickness of corrosion products in this model) is found 
by substituting Eq. (2-8) into Eq. (2-6): 


n 


c 


r 




(2-9) 


Thus, as the temperature increases, both t c and n c /A decrease. 

This is reasonable since less catalytic material would be required to obtain 
a given gas generation rate as the temperature increases. The value of 
n c -v- 6 x 10"® 16-mole ( n c/A ^ 3.6 x 10"^ lb-mole/in^) which was observed 
as the approximate value of the critical gas content in the temperature 
range 160-1 79°C agrees reasonably well with Eq. (2-9). At 170°C, Eq. (2-9) 
yields n £ /A * 4.3 x 10"^ lb-mole/in 2 . 



■ — * 


13111-6060-RU-00 


2.3 Conclusions and Recommendations: 

The behavior of the gas evolution in methanol/stainless-steel 
heat pipes was found to depend on the amount of gas per unit area generated 
during accelerated testing. Below a critical value n c /A given by Eq. (2-9), 
or critical time t_ given by Eq. (2-8), the time dependence is explained best 

v 

by a parabolic function, indicative of the growth of a passivating film of 
corrosion products. The data in this region can be correlated with a model 
of film growth resulting in Eq. (2-6), which contains parameters to be 
determined by experiment. Above the critical point given by Eqs. (2-8) 
and (2-9), a linear time dependence was observed. In this region the 
data can be correlated with a model of catalytic decomposition of methanol 
to formaldehyde and hydrogen on the surface of the corrosion products, 
resulting in Eq. (2-7). No flow rate dependence was found within the 
accuracy of the data. In application to other types of methanol /stainless- 
steel heat pipes, the gas generation may vary depending on the type of 
stainless steel, the purity of the methanol, and other factors; but it is 
expected that the behavior could be explained by the same form of the Eqs. 
(2-6) and (2-7), with associated critical values (2-8) and (2-9), only 
the value of the parameters may change. 

Based on the results of this and the previous study [2-1] it 
appears that this method of accelerated life testings has a broad 
applicability to heat pipe systems, even when not a great deal is 
known concerning the actual gas evolution mechanisms. This method of 
accelerated life testing can now be applied to other important types 
of heat pipes with good probability of success. 
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3.0 6ASPIPE EXTENSIONS 

The TRW GASPIPE computer program, which was developed on the 
present contract and documented in Ref. [3-1], "User's Manual for the TRW 
Gaspipe Program", has been extended to include two condenser sections. 

Since the new version of the program is described *.n detail in the rented 
user's manual, (Ref. [3-2]), it will suffice here to summarize the 'irograms 
capabil ities. 

The motivation for extending the original GASPIPE pro ram 
is the frequent design of heat pipes with a second condenser section r*vi,g 
as low an axial conductivity and effective sink temperature as possible. 

Such sections are used adjacent to the gas-reservoir entrance to minirize 
the partial pressure of the reservoir vapor, thus allowing smaller 
reservoirs and/or a tighter control band. In the case of a wicked cola- 
reservoir design, the low-conductivity section minimizes heating of the 
reservoir by axial conduction from the active portion of the condenser. 

In the case of a non-wicked hot-reservoir design, wherein the partial 
pressure of the reservoir vapor is set by the temperature of the nearest 
portion of saturated wick, the low-conductivity section minimizes warming 
of this portion. 

The original GASPIPE program cannot treat such situations if 
the sink temperatures of the two condensers are unequal and the gas front 
develops in the condenser farthest from the reservoir. In addition, the 
original program is limited to either two condensers or one condenser and 
one adiabatic section, and many applications require two condensers and 
an adiabatic section. The extension of the program remedies these 
deficiencies. 

The situation to which the revised program, GASPIPE 2, is applic- 
able is depicted schematically in Figure 3-1. In this figure, the axial 
conductivity of condenser 1 is minimized by the use of slotted fins. All 
parameters including the sink temperature that describe condenser 1, can 
be set differently from those describing condenser 2. Thus the fins attached 
to condenser 1 may have a special coating to achieve a low effective sink 
temperature, or the heat-pipe wall itself may be thinner or even be a different 
material to achieve a low axial conductivity. A typical condenser cross- 
section is shown in Figure 3-2. The program Is not limited to the geometry 
shown, non-circular and non-axi symmetric configurations can also be studied. 
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In addition to the extensions, the new GASPIPE 2 program retains 
all the capabilities of the original GASPIPE program, which allow one to: 

o Calculate the wall-temperature profile along a gas 
loaded heat pipe. 

o Calculate the amount of gas loading necessary to 
obtain a desired evaporator temperature at a 
desired heat load. 

o Calculate the heat load versus the evaporator 
temperature for a fixed amount of gas in the pipe. 

o Calculate the heat and mass transfer along the 
pipe, including the vapor-gas front region. 

o Calculate the heat leak when the condenser is 
filled with gas. 

o Calculate whether or not freezing occurs in the 
condenser and, if so, at what rate. 

o Determine the information required to size the gas 
reservoir of gas-controlled heat pipes. 

The program contains numerous reservoir options that allow it to be used 
for hot or cold passive-control as well as heated-reservoir active- 
control heat pipes. 

GASPIPE 2 Is based on a one-dimensional steady-state analysis 
that results in two simultaneous first-order ordinary differential 
equations that govern (i) the mole fraction of the noncondensable gas 
and (ii) the vapor velocity. These are Integrated with the forth-order 
Runge-Kutta routine. The addition of the second condenser section with a 
different sink temperature required extensive modification to the original 
program. Previously, a single solution curve was obtained by numerical 
integration that could oe translated within the condenser until a position 
corresponding to the required heat flow or the required inventory was 
achieved. With a step chaiige in condenser properties that occurs with 
two condensers, translation is no longer possible. An Initial solution 
Is obtained for the case of complete gas blockage of the heat pipes by 
first Integrating once backward from the adiabatic section with properties 
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of condenser 2 and then repeatedly forward with properties of condenser 1 
until proper conditions are met at the intersection of the two curves. The 
two curves are then translated until their intersection coincides with the 
boundary between condensers 1 and 2. With the initial gas-blocked profile as a 
base, a solution curve for a non-gas-blocked region is obtained by select- 
ing a point along the gas-blocked curve, slightly perturbing the value of 
the solution at that point, and then integrating towaru the evaporator. 

The solution curve thus generated automatically see*s a condition of no 
gas blockage. Solutions are repeatedly obtained in this way, but with 
different starting points along the gas-blocked curve, until either the 
calculated amount of noncondensable gas or the total heat rejected agrees 
with the specified amount. In some cases, the correct solution is obtained 
by repeated integrations from the beginning condenser 1 with suitable adjustments 
of the initial conditions. For details of the program, refer to the user's 
manual already cited. 
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4.0 PRIMING STUDIES 

Many heat-pipe applications require both the high performance 
of an arterial design and the variable-conductance characteristics obtained 
through the use of noncondensable gas. The presence of gas, however, can 
interfere with the priming and operation of arteries. Part of our research 
task was to study arterial performance in the presence of gas. 

During arterial priming, noncondensable gas can result in 
bubble entrapment. At TRW, we have developed an approach to vent gas 
through holes in a foil-walled section of the artery at the evaporator end. 
Our research task calls for detailed study of this new solution of arterial 
priming with a glass heat pipe. As will be described presently, such an 
apparatus was fabricated that closely approximates actual heat-pipe 
configurations currently in use. 

Noncondensable gas also causes problems in the operation of 
arterial heat pipes with ammonia. Pressure fluctuations which only appear 
when there is some condenser blockage by nor condensable gas, result in 
depriming cf arteries. Our research task also calls for both a theoret- 
ical and experimental investigation of the mechanism of these fluctuations. 

4.1 Theoretical Study of Pressure Fluctuations in a Gas - 

Controlled Ammonia Heat Pipe 

Pressure fluctuations in gas-controlled ammonia heat pipes 
were first reported by Edelsteln, Roukis and Loose (Ref. [4-1]). They 
observed 0.12 to 0.14 psi irregular fluctuations with a period ranging 
between 0.25 and 0.50 minutes in an ammonia heat pipe containing nitrogen 
gas for control of the active condenser length. The pipe was approximately 
one inch in diameter with a 46 inch long evaporator and an equally long 
condenser. The magnitude of the fluctuations increased somewhat with 
heat load, but was independent of tilt Their magnitude was sufficient 
to deprime the tunnel artery, which was the basis for their high perform- 
ance (150,000 watt-inch) heat pipe. With no control gas present, they 
reported a marked decrease in amplitude and "number" (frequency?) of 
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oscillations so that the pipe performed as desired. 

Similar problems with the priming of a high-perfomance ammonia 
heat pipe were encountered by TRW (Ref. [4-2]). A pressure transducer was 
installed and indicated excessive pressure fluctuations. In this case 
the fluctuations, while still rather irregular, were much more nearly 
periodic in nature with a period of approximately 1.1 to 1.2 minutes 
and an amplitude of approximately 0.31 Ib^/in^ at low heat load. Data 
were obtained with and without additional mass added to the evaporator 


as follows: 

Mass of Aluminum 
on Evaporator 

lb_ 

m 

Control Gas 

Steady Power 
Watts 

Pressure-Fl uctuati on 
Amplitude 

lb f /in 2 

0 

None 

100 

< 0.02 


Argon 

100 

> 

o 

"o 

--4 


1) 

3 

-0.45 

3.5 

II 

100 

< 0.07 


II 

S 

< 0.33 

28.0 

II 

100 

< 0.15 


II 

3 

< 0.010 


Contrary to the experience reported by Edel stein, the most 
severe oscillations occurred at low heat loads. However, in these tests 
the condenser sink was massive, and the severe oscillations occurred when, 
at low heat loads, the gas front penetrated Into the adiabatic section 
whose mass was only that of a 0. 50-Inch outside diameter tube with 0.028- 
inch thickness and the diametral wick. 


4.1.1 Possible Oscillation Mechanisms 

The long period of the oscillations, 15 to 70 seconds, strongly 
suggests a mechanism limited by heat or mass-diffusion rates. The 
possibility of gravity waves playing a part, for example, may be ruled 
out by considering the time scale for such a mechanism 


t , 36 

gravity 



4 0.1 Sec. 
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where p is the density of the gas, D Is the pipe diameter, and g is the 
gravitational acceleration. Such considerations would rule out with 
even greater force acoustlc-resonar.ce phenomena whose period would be 
on the order of the pipe length divided by the acoustic velocity. In 
comparison, consider typical characteristic times for diffusion 
phenomena shown in Table 4-1. We will see presently that the theoretical 
period of the oscillation is of the same order of magnitude (within a 
factor of 4) as a characteristic delay time. Therefore, the times most 
in agreement with the experimental ovservatlons are those for phenomena 
1 , 2, 3 and 5 in Table 4-1. 


TABLE 4-1 

CHARACTERISTIC TIMES FOR 
DIFFUSION PHENOMENA 



PHENOMENON 

FORMULA FOR 
CHARACTERISTIC TIME 

NUMERICAL VALUE OF 
CHARACTERISTIC TIME 

1 . 

Mass diffusion in the 

(pipe diameter) 2 

10 Sec. 


gas phase through one 

(diffusion coef.) 



inside diameter 



2. 

Mass diffusion in the 

(gas-front length) 2 

100 Sec. 


gas phase through the 
gas front 

(diffusion coef.) 


3. 

Mass diffusion in the 

(wick thickness) 2 

50 Sec. 


liquid phase through 
the wick thickness 

(diffusion coef.) 


4. 

Heat diffusion in the 

2 

(groove depth) 

0.1 Sec. 


liquid phase through 
groove depth 

(thermal dlffuslvity) 


5. 

Heat diffusion in the 

(wick thickness) 2 

10 Sec. 


liquid phase through 

(thermal dlffuslvity) 



the diametral wick 
half thickness 
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Three poss^le disturbances can be conceived on the basis of 

this table: 

A. A disturbance based upon the time lag between a convective 
distortion of the gas front and Its re-establishment by 
diffusion. 

B. A disturbance based upon a convective distortion of the 
gas front and the consequent effect upon dissolved gas in 
the condensate which in turn triggers nucleation in the 
evaporator. 

C. A disturbance based upon the time lag between the filling 
of a void in the diametral wick with liquid and the 
formation of a bubble by boiling in the wick where it 
contacts the wall. 

These three correspond with Items (1 or 2), (3), and (5) 
respectively in Table 4-1. Item ( 4 ), which corresponds to periodic bo 
boiling in the grooves. Is seen to have too fast a characteristic time. 

Further consideration of the time necessary for the liquid to 
flow from condenser to evaporator would seem to eliminate disturbance B 
above. Further, while a diametral wick was used in the TRW experiments, 
a spiral arter> was used in work of Ref. [4-1]. The feed wicks of that 
pipe do not seem to have a configuration favoring foiling within them, 
thus disturbance C seems unlikely. Disturbance A above is thought most 
likely, and we now consider a simplified model of It. 

4.1.2 A Simplified Model 

In a heat pipe a pressure change dP is caused by (or accompanied 
by) a temperature change dT given by the Clausius-Clapeyron relation, 
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where h^ and M respectively are the latent heat and molecular weight 
of the working fluid, and^?is the universal gas constant. The observed 

_3 

fluctuates in dP/P are on the order of 10 ; hence, the fluctuation in 

-4 

dT/T is on the order of 10 . The pressure and temperature fluctuation 

are also accompanied by a prompt movement dz of the gas front. If A c 
is the cross-sectional vapor area of the condenser, and V is the total 
gas-blocked voluvie that contains N moles of vapor and gas at temperature 
T $ , then 


NflT. 


dz A„ ■ -dV --d( — p—2.) = V ^ 
or, using (4-1), we write 

-iSf f • 


(4-2) 


c ' ' s 

The movement of the gas front by dz opens up new heat-transfer area so 
that the heat-transfer rate Q lost in the condenser goes up by dQ-j: 


dQ 1 * UPdz (T - T $ ), (4-3) 

or, using (2) , we write 

dQ, - UP (T-T,) ^ f • (4-3) 

Here, U is the heat- transfer coefficient and P is the condenser 
perimeter. There is also a slightly augmented heat transfer due to 
the increase in T - T , 

dQ 2 * UPz dT (4-4) 

where z is the active condenser length. We imagine that the latter is 
prompt in terms of a few seconds while the former may be delayed by a 
time x of some tens of seconds or a minute due to the slowness with 
which the gas diffuses out of the way of the condensing vapors. 

The heat loss from the condenser In time dT Is O c ■ dQ^ + dQ 2 , 
where we can write Q c as 
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■ < T • V ^ 



dR o 

ar 


(T(t-r) - T q )3 . 


(4-5) 


Here, T Q is the mean vapor temperature, and from Eqs. (4-2), (4-3) and 
(4-4), we have (to fir.»t order in T - T ) 



( 4 - 6 ) 


( 4 - 7 ) 


In order for the evaporator temperature T to rise by dT a quantity 
of heat dQ is 'romptly needed to heat the evaporator, adiabatic, and 

caC 

condenser masses: 


dl W = f m e c e + ra a c a + m c (z/ ’ L c) c cJ dT 

where m and c are the mass and specific heat of the evaporator, adiabatic 
or condenser sections, which ar° distinguished by the subscripts ( ) e , 

( ), and ( 1 respecLi .aly , and is the total condenser lenqth. There 
is also a delayed need of heat dQ z needed to warm the segment of condenser 
brought into service: 


dQ z ■ (m c /L c )c c (T - T $ ) dz. 


or, using (4-2), 

W^c>c c <T - V j(- f . 

The steady electrical heating of the evaporator goes into supplying this 
heat as well as the heat losses [Eq. (4-5)3 during a temperature fluctuation], 
that is 
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+ (T - V r n + 1 
0 0 


dR ft 

(T(t - t) - T 0 )]dt. 


Denote 


T p ■ R 0 [ Ve + V. * "c ( V L c )c c 3 - 

T - T c „ h- M 
T <J = R o [(m c /L c lc c ( -^ — 1 


T c - T s * R o «• 


B = (T - TJ 


1 dR o 


s' r ar - 


T ‘V 


r . t - t 0 . 


The heat balance equation becomes 


(4-8) 

(4-9) 

(4-10) 

4-11) 

(4-12) 


0 = 


dT’ 

ar 


T dT ’ 

T d ar 


+ T* (t) + BT'(t - t) 


+ B I Mllll. (4-13) 

0 s 

Eq. (4-13) is a differential -difference equation which are discussed 
in Refs. [4-3] and [4-4]. It einbodies a prompt effect of the thermal 
capacities of the evaporator adiabatic and condenser sections in resisting 
temperature change (the first term), a delayed effect of the thermal capacity 
of the condenser exposed when the front moves (the second term) , a prompt 
effect of Increased heat losses in the working portion of the condenser when 
the vapor temperature rises (the third term), a delayed effect of Increased 
heat loss due to movement of the front In exposing more condenser surface (the 
fourth and fifth terms). This last factor is represented by the first- 
order linear term BT' (t - x) representing exposure of new condenser 
surface to the old vapor temperature T Q and a second-order, nonlinear 
term [B/(T Q - T $ )] T(t) T* (t - x) representing the Increased heat 
transfer from new surface due to the rise T' in vapor temperature above 
T q . This latter Is clearly small and will be neglected when the equation 
is examined for what it says about the period of the disturbance. 
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In order to say anything about the amplitude of the disturbance, 
it is necessary to retain nonlinear terms. That the only nonlinear 
term shown in Eq (4-13) is the dominant one is not clear. In the course 
of deriving Eq. (4-13) many inherently nonlinear terms were neglected; 
for example, the change in gas volume V due to movement of the front or 
to the change in temperature T (taken to be T ) in the equations for 
the B and terms [Eqs. (4-7) and (4-11) and (4-9)]. Since the factor 
of concern is V/T and V decreases upon increase in T, the two effects 
are additive rather than self-cancelling. Furthermore, in the tests of 
Ref. (4-1), spray cooling was used, and radiation is often an important 
heat-transfer mechanism in rejecting heat from the condenser in actual 
applications. Where either evaporation (spray cooling) or boiling or 
radiation acts at the condenser, an increase in T - T gives more than 
a proportional increase in the heat rejected, and hence another source 
of nonlinearity arises. A mitigating factor in favor of retaining the 
nonlinear term of (4-13) for an amplitude analysis is that it has much 
the same form of a "stiffening term" of a spring-mass system. 

4.1.3 Analysis of the Equation 

Much information can be obtained from Eq. (4-13) if we simplify 
it by neglecting the nonlinear term and considering the special case 
t d = 0. In this case, we have 

T p 3T t + T ' (t) + B T '( t “ T ) a °* (4-14) 

We seek a solution of the form 

T' = A e* t , 
where c = a + i a.. 

If a is positive, our linearized theory predicts that the solution will 
become infinite. Actually, the fluctuations will be limited by nonlinearities 
or other small effects that we have neglected. If o is negative, the 
solution is damped and we do not expect fluctuations. 
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The characteristic equation is obtained by substituting 
Eq. (4-15) into (4-14): 




Ct 


- o. 


(4-16) 


or in terms of its real and imaginary parts, 


T p a + 1 + B e" aT cos air = 0, 


(4-17) 


oj T p - B e' aT sin u>t = 0 


(4-18) 


We first ask if there are any non-oscillatory solutions U = 0 ) 
that are unstable (a > 0). By setting «o ■ 0 in (4-17), we obtain 


a = -1 -Be 


-OT 


(4-19) 



Since B > 0, we deduce that if the solution is non-oscillatory, then it 
is stable. 

We now look for an unstable (a > 0) oscillatory solution 
(w > 0). From (4-17) we obtain 


u)T = COS 


■1 /-I -TjM 

\ B e" aT / ‘ 


(4-20) 


According to this equation, fluctuations originate (u = 0) at 
B = 1, in which case they have a period 2w/u = 2t. This result applies 
to a hypothetical heat pipe with zero mass since we must have r p = 0 to 
satisfy Eq. (4-18) and hence allow such a solution. In the opposite 
limit of a massive heat-pipe (t p » l), oscillations originate (a • 0) 
only ’f B » 1, in which case we see from Eq. (4-20) they have a period 
2ir/w ■ 4 T> anC | f rom Eq. (4-18) the necessary value of B is 
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Since B is proportional to the gas-blocked volume, the fact 
that B has a critical value at the onset of unstable oscillations suggests 
that a critical gas-blocked volume V cr exists. For a zero-mass heat pipe 
(tp = 0), which has a critical value B * 1 , we have from Eqs. (4-11 and 4-12) 



A c 


S M 



(massless heat pipe). 


(4-21) 


whereas for a massive heat pipe (t >> 1), which has a critical value 

p 

B = u)Tp, we have from Eqs. (4-11), (4-7) and (4-18) 

(massive heat pipe) 


Fluctuations in a heat pipe with gas-blocked volume smaller than its 

critical volume are damped, whereas the converse is true if the gas- 

blocked volume is larger. From the above expressions we see that a 

heat pipe becomes increasingly more unstable as the active condenser 

length z 0 approaches zero. In fact, our simplified theory predicts that 

the massless heat pipe will be unstable even with a small gas reservoir since 

as z Q -*■ 0, V cr -*-0. A low sink temperature is also a destablizing factor, that 

is, V cr is more likely to fall below the gas-blocked volume if T s << T 0 . 

For a massive heat pipe, the critical volume is inversely proportional 
to the time lag. Therefore, a short time lag favors stability. The 
time lag can be decreased by increasing the gas diffusivity if either 
Item 1 or 2 of Table 4-1 governs or by decreasing the wall thickness or 
conductivity if Item 2 governs. 


When is non-zero in Eq. (4-13) much the same reasoning 
applies, but V cr is increased. Approximately, the critical volime is 
given by 



* <T P * Td ) T s 

2 < Wg" ) V T S 


(4-23) 
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Since both T p and are inversely proportional to z 0 [Lqs . (4-6), (4-8) 
ana (4-9)], V cr remains finite as z Q approaches zero. The effect of 
including the ^ term in the linearized analysis does not change the 
preceeding qualitative discussions; it merely increases the size of the 
critical volume. 

The theoretical model may be summarized as follows: 

o A diffusion time-lag phenomenon is thought to govern. 

The re-establishment by diffusion of a convectively 
distorted vapor-gas front is thought to be the mechanism 
of the lag. 

o The period of instability ranges from 2 to 4 times 
the diffusion time lag. 

o A gas-blocked volume greater than a critical volume 
V cr necessar y f° r fluctuations to occur. V cr 
increases with increasing thermal mass and decreases 
with decreasing sink temperature. 

o The heat pipe tends to be more unstable at low 
power levels when the active condenser length is 
small. A massless heat pipe is always unstable at 
low power, but thermal mass can stabilize at all 
powers. 

o A high mass-diffusivity control gas such as helium 
and low axially conducting condenser wall such as 
stainless steel favor stability. 
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4.2 Experimental Measurements of Pressure Fluctuations 

The theoretical model of pressure fluctuations is based on a 
diffusive time lag in the establishment of a vapor-gas front in response 
to a change in vapor temperature. As a crucial qualitative test of the 
model, a series of pressure-fluctuation measurements were made with 
a heat pipe containing as the control gas first argon, which has a 
relatively low diffusivity and then helium, which has a relatively high 
diffusivity. The diffusive time lag for helium should be much shorter 
than the lag for argon; thus if the postulated model is correct, the 
change of control gas should have a marked effect on the fluctuations. 

The specifications of the test setup are contained in Table 4-2. 
The heat pipe used 'or the experiment was a stainless-steel two-artery 
slab-wick design, 62.75 inches long and 0.5 inches in diameter. The arteries 
were never primed during the tests, but sufficient heat-transfer capacity 
in excess of 80 watts was obtained in a level orientation. 

Pressure fluctuations were measured with a 5-psia 
differential pressure transducer, both sides of which were connected to 
the feed tube in the heat-pipe reservoir. One leg of the connection is 
fitted with a valve that is closed when measurements are being made and 
open otherwise. T he output from the transducer is recorded with a strip- 
chart recorder. Runs were made at low power, in which case the gas front 
extends into the adiabatic section, and at intermediate power, in which 
case the front is in the condenser. 

Typical fluctuation measurements are shown in Figs. 4-1 and 

4-2. With argon at the low-power of 3 watts (Fig. 4-1), the pressure 

trace varied slowly and randomly, typically at a rate of - .05 psi /minute 

except at discrete points when the trace suddenly recorded a marked rise 

or fall of approximately .35 psi in less than half a minute. Such an 

« 

isolated surge is shown in the Figure (4-1). Only two surges were 
recorded in a 31 minute interval. At a higher power of 10 watts, the 
front was still in the adiabatic section, and while the trace is similar 
to that at 3 watts, six isolated surges were recorded in 44 minutes. These 
surges alternately Increased and decreased. In comparison, at the low 
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Table 4-2. 

Specifications of the Test Setup 
for Pressure Fluctuation Measurements 

Keat-pipe dimensions (all materials stainless steel): 

Outside diameter 
Wall thickness 
Internal grooves 
Evaporator length 
Adiabatic length 
Condenser length 
Cross-sectional area 
of vapor space 
Reservoir volume 
Arteries (2) 

Wick 

Test conditions: 

Thermal resistance between 
condenser saddle and heat sink - 1/32 x 1.5 x 25 inch Teflon. 

Condenser saddle - .1275 lb^inch of aluminum. 

Set point - with the sink and reservoir 

at 10°F, sufficient gas was 
added to give a 70°F vapor 
temperature at 67 watts. 


0.5 inches 
0.028 inches 
150/inch 
23 inches 
16.75 inches 
25 inches 

.1064 square inches 
8.72 cubic inches 
.063 inch I.D. 

.050 inch slab 





41 




FIGURE 4-1. Pressure fluctuations at low power (3 watts) with the vapor-gas front in the adiabatic region 
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FIGURE 4-2. Pressure fluctuations at higher power (67 watts) with the 
vapor-gas front In the condenser region. 
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power of 3 watts with helium as the control gas, no surges were observed. 
However, as seen in Fig. 4-1, the trace is characterized by fairly regular 
oscillations of about 7 cycles per minute around a randomly varying mean 
pressure. 

At the higher power, 67 watts, in which case approximately one 
fourth of the condenser is active, the situation is somewhat reversed. 

As seen in Fig. 4-2, no surges were observed with argon, however, some 
were observed with helium. Also, the oscillations with helium are more 
pronounced and of a higher frequency (17 cycles/minute) at the higher 
power. 

The primary conclusion of the experimental measurements is 
that the type of control gas, as the theory predicts, has a marked effect 
on the fluctuations. This strongly suggests that a diffusion time lag plays 
a crucial role in the fluctuation mechanism. We cannot conclude, however, 
that pressure fluctuations are well understood. For example, we observed 
pressure surges, that have no counterpart in the theoretical model. Further 
experimental research is required to verify the key physical parameters that 
effect the fluctuations and to assess the magnitude of fluctuations that 
can be tolerated without Interfering with arterial operation. It may then 
be possible to design an arterial variable-conductance heat pipe with 
ammonia as tne working fluid. 

4.3 Priming Studies with a Glass Heat Pipe 

A method has been developed at TRW to prevent entrapment of 
a gas bubble during arterial priming. The noncondensable gas is vented 
through capillary-size holes in a foil-walled section of artery at the 
evaporator end. Liquid cannot plug the holes in the foil and thus prevent 
venting, because the foil is sized so thin that the menisci on either side 
of a potential plug would coalesce. 

Part of our research task was to fabricate a glass heat pipe 
and use it to study priming with this new approach. Several aspects of 
arterial priming can effectively be studied with a glass heat pipe. For 
example, venting gas at the evaporator end of the artery relies on priming 
proceeding uniformly from the condenser end to the evaporator end. If 
a bubble is trapped at some position along the artery other than at the 
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evaporator end, then it would first have to be convected to that end 
before it could vent. Convection of a bubble for various heat loads 
can be observed directly with the glass heat pipe. Another process to 
be studied is priming under a load. For a given evaporator elevation, 
we can use the glass heat pipe to find the maximum heat load under which 
priming will occur. Other details of the priming process can be studied 
such as the minimum capillary hole size in the foil that will still vent 
gas, and whether in some circumstances, the foil -walled portion of the 
artery can be flooded by a pool of excess liquid and hence be unable to 
vent gas. We were unable to answer these crucial questions because funds 

on the contract were exhausted before the studies could be carried out. 

Fabrication of the apparatus is essentially complete. As shown 
in Fig. 4-3, the main structure consists of a 0.752 in. O.D. and 0.515 ir.l n. 
glass tube, 43 inches long, with machined stainless-steel end fittings. 

The heat pipe is held in a plexiglass cradle, which also serves as a 
safety shield when ammonia is used as the working fluid. 

The heat pipe is designed for study of the arterial/slab-wick 
configuration. Cross sections in the evaporator and condenser sections 
are shown in Fig. 4-4. The artery is mounted on one side of the slab 
wick, and a 1/8-inch diameter stainless-steel sheathed heater with a 12-inch 
heated section is inserted into a double-layer screen casing that is spot 
welded to the other side of the wick. The unheated end of the heater passes 
through and is brazed into the evaporator end cap. In the condenser region, 
a stainless-steel cooling loop extends along the opposite side of the wick 
as the artery. It passes through and is brazed into the condenser end cap. 

A 1/16-inch O.D. stainless-steel sheathed thermocouple extends into the 
adiabatic section so the vapor temperature can be monitored. A feature of 
this design is that it closely models actual arterial /slab-wick heat pipes 
such as the ones to be used in the rocket-flight experiment as described in 
the next section. 
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5-0 SOUNDING-ROCKET EXPERIMENT 

The GSFC sounding-rocket heat-pipe experiment scheduled for 
1974 provides a unique first opportunity for a zero-grav'ty test of 
the mechanism of menisci coalescence to vent noncondensable gas during 
arterial priming. Part of our research task was to design, fabricate 
and test two research heat pipes (and two spares) for the flight 
experiment. 

The experiment is especially challenging in that only six 
minutes in zero gravity are available for priming and then applying 
a heat load to verify that a primed state was achieved. In addition, 
there are uncertainties as to whetner the so-called "pogo effort" uni 
interfere witn priminq. If, at the instant of rocket-engine shut-down, 
a surge of excess liquid arrives at the evaporator before the artery 
primes, the evaporator end of the artery may fill with liquid and prevent 
venting. Our task is not only to test whether or not the heat pipe primes, 

but also to obtain some crucial information on the internal conditions of 
the heat pipe during the experiment. Thus, each heat pipe has one of two 
otherwise identical arteries instrumented at the evaporator end with a 
miniature thermistor. The thermistor, which is electrically heated, is 
markedly cooler when submerged in liquid than when in vapor. As such, it 
detects the presence or absence of liquid in the priming foil. In this 
way, the research heat pipes are designed to overcome a nagging difficulty 
of arterial heat pipes, that of diagnosing the cause of unsuccessful priming. 
Priming failure can have several causes, however, they generally fall into 
one of two categories listed below. 

Two Categories of Priming Failure 

1. Evaporator stress too high for priming due to 

o insufficient fluid in the heat pipe 
o too large a residual heat load during priming 

o too large of a hydrostatic load due to adverse 
orientation in an acceleration field 
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2. Entrapment of noncondensable gas due to 

o failure to vent the gas through the priming foil 
o a gas bubble at some location in artery 
ot.ier than the evaporator end. 

If, during the priming period, liquid no£ detected at the 
evaporator end of the artery, then priming failure is attributed to the 
first category, excessive eva r Drator stress.^ If, on the other hand, 
liquid detected, we conclude that the stress was low enough for 
priming to take place, and a subsequent failure of the uninstrumented 
functional artery is attributed to the second category, entrapment of 
noncondensable gas. 

The thermistor instrumentation also provides a means of 
verifying that a sufficient heat load is applied when we test for the 
primed state of the functional artery. If the thermistor detects that 
liquid recedes (which it should, due to the large pore in the 
instrumented artery), then the heat load is sufficient to empty the 
functional artery if pricing is unsuccessful . If subsequently a 
burnout does not occur, then the primed state of the functional artery 
is verified. 

The diagnostic logic is summarized in Figure 5-1. The heavy 
line represents successful priming and verification of such. 

5.1 Asc r iption of the Heat Pipes 

As shown in Figure 5-2, the heat pipes are of an arterial/slab- 
wick configuration. Their specifications are listed in Table 5-j. The 
two heat pipes to be used in the flight experiment differ only in the 
details of their priming foils. For the heat pipe designated X-2, the 
active length of the priming foil is 1.5 inches, the holes are in spiral 
rows, and the end on the functional artery Is crimped closed, whereas 
for the heat pipe designated X-l , the active length Is 0.375 inches, 

1. A relative large pore of comparable size to the artery diameter is 
provided just behind the thermistor to Insure that gas vents from 
the Instrumented artery. 
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Thermistor detect NO 

liquid in end of artery 

during priming period? 


Heat-pipe burn-out 
when load applied? 


Instrumentation 

malfunction. 


Failure due to 
excessive evaporator 
stress during 
priming period. 


Thermistor detect 
liquid leaving 
end of artery 
when load is applied? 


Heat-pipe 

burn-out? 


Heat-pipe 

burn-out? 


Heat load may not 
have been high 
enough for 
definitive test. 


Instrunentation 

malfunction. 


Verification of 
primed state of 
functional artery. 


Failure due to 
bubble entrapment. 


FIGURE 5-1. Sumary of diagnostic logic for the flight experiment. 















FIGURE 5-2. Configuration of the flight heat pipe. 
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TABLE 5-1 

SPECIFICATION OF THE RESEARCH HEAT PIPES 

Tube (304 stainless-steel): 

Length - 36 inches 
Outside diameter - 0.5 inches 
Wall thickness - 0.028 inches 
Internal threads - 150 per inch 

Slab-wick (304 stainless-steel X-13 felt metal): 
Thickness - 0.050 inches 
Fiber diameter - 0.00085 inches 
Porosity - 84% 

Arteries (304 stainless-steel 150-mesh screen): 

Inside diameter - 0.063 inches 

Height from bottom of slab wick - 0.170 inches 

Priming foils (304 stainless-steel): 

Thickness - 0.0005 inches 
Diameter of holes - 0.010 inches 
Hole spacing ir. row - 0.032 inches 
Row spacing - 0.075 inches 

1 1.50 inches for Configuration X-2 
0.375 inches for Configuration X-l 

! 17° spirals for Configuration X-2 
straight for Configuration X-l 
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TABLE 5.1 (cont.) 

Thermistor (Veco, part No. 32A7): 

Bead diameter - .013 inches 
Leads - .001 inch platinum 
Resistance at 25°C - 2000 ft 

Flight heat sink: aluminum block 12" X 1.32" X 1.32" 
Flight heater: 10 inches long and beginning 2 inches 
from the evaporator end 

Flight thermistor locations: 

Thermistor No. 


4 

5 

6 

7 

8 
9 


Distance from evaporator end 
3 inches 
7 

16 

20 

25 

30 

35 

on heat sink 
on heat sink 
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the holes are in straight rows and the end is closed with a plug. The 
X-2 configuration priming foil is a design similar to that used on the 
Canadian Technology Satellite (CTS) program. The X-l configuration 
priming foil is an experimental design that represents an attempt to 
improve the ruggedness and minimize the chance of trapping a gas bubble 
in the last instant of priming the foil itself. 

The instrumentation consists of a bead- type thermistor that 
is held at the apex of the conical end of a cylindrical teflon plug, 
(.073-inch diameter, .4-inch long) which in turn fits into the end of 
one of the priming fGils. The .001 -inch platinum leads from the bead 
pass through two 0.008" diameter holes that run the length of the pl'jg. 

The plug is held in the end of the priming foil by a band of heat-shrink 
tubing. The thermistor leads are connected to the prongs of an electrical 
feed through in the heat-pipe end cap with 0.007-inch nickel wire. 
Connections are made with spot welds. The nickel wire is mechanically 
secured to the outer end of the teflon plug with another band of heat- 
shrink tubing. 

5.2 Fabrication Highlights 

Some problems arose during fabrication, which are discussed 
here so they can be avoided in future projects. 

We originally planned to use a miniature hot-wire probe 
instead of a thermistor as the liquid sensor. During the breadboard 
test, however, we discovered that methanol attacked the solder that 
secures the hot-wire element to the support prongs. After a 12-hour 
submersion in methanol, the hot wire detached. At this point we decided 
to use a thermistor. 

Because of the incompatibility of methanol and solder, we 
could not solder the electrical feed-through into the heat-pipe end cap. 

We tried heli-arc welding, but the heat generated resulted in the cracking 
of the glass insulation. Laser welding, which generates little heat, was 
also tried, however, we found that the weld bead had microscopic cracks 
that resulted from the dissimilarity of inconel of the feed-through and 
the stainless-steel end cap. The solution was the use of a pure-lead 
braze after both parts were first gold plated. As a final assembly step, 
the end cap was laser welded to the heat pipe. 
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Another problem occurred as the result of vacuum firing. Since 
we did not want to subject the thermistor installation to the vacuum- 
firing temperature, we planned to vacuum fire before the thermistor was 
installed. The evaporator end of the wick was left extended approximately 
one inch from the tube, the heat pipes were vacuum fired, and the thermistor 
was installed. We then planned to pull the wick the final distance into 
the tube in the usual manner, however, vacuum firing resulted in seizure. 
Only after pulling with excessive force, were we able to pull a wick into 
one of the heat pipes. The other three heat pipes were completed by 
cutting off the condenser end of tube and welding a short piece of tube 
to the evaporator end with the use of a sleeve over the joint. 

5.3 Fill Determi nation 

An initial fluid fill of .0390 lb, of methanol was calculated 

m 

by TRW's MULTIWICK computer program to give the maximum heat load under 
which priming will occur in zero gravity without a vapor-space slug. 
MULTIWICK also calculated that .0390 lb, of methanol would be more than 
sufficient to prime in earth gravity. In the computer calculations, 
however, we did not take into account an effect that was observed in 
a sub-scale glass-tube priming test; that the stress required to prime the 
foil -walled portion of an artery is somewhat lower than that required to 
prime the screen-walled portion. In initial tests, priming was difficult 
to achieve. The problem was solved by increasing the amount of methanol 
to .0487 lb m> This amount corresponds to the maximum in zero gravity 
that will not result in a vapor-space slug with zero heat load. Under 
load, however, some slugging will occur. 

Since the experiment concerns priming in the presence of non- 

_g 

condensable gas, 1.87 x 10 lb moles of a mixture of 90% nitrogen and 
10% helium were added to each heat pipe. This amount results . 
approximately one-inch of condenser blockage. 

5.4 Steady-State Capacity Tests 

Steady-state capacity tests were carried out and the results 
are summarized in Figure 5-3. The theoretical predictions based on 
MULTIWICK are shown as well. Each of the four heat pipes were tested 
with primed arteries at 1.0-inch and 0.25-inch elevations, however. 
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only one was tested with failed arteries. Agreement between the theoretical 
predictions and the experimental results are generally good except in the 
case of 0.25-inch elevation where predictions are generally too high by 20%. 

5.5 Priming Tests with the Thermistor Instrumentation 

The electronic circuit for the thermistor is particularly 
simple; it consists of a 3000-ohm resistor and a 28-volt power supply in 
series with the thermistor. The circuit output is the voltage across the 
thermistor. The output voltage is approximately 3.7 volts when the arteries 
are primed and 1.2 volts when deprimed. Sometimes, however, in the primed 
state boiling occurs at the bead, in which case the output oscillates 
rapidly about some mean voltage greater than 3.5 volts. In the deprimed 
state, discrete intermediate output voltages occur that are attributed to 
partial cooling of the thermistor bead by small liquid fillets around it. 
With sufficient evaporator stress these fillets are not replenished from 
neighboring liquid and they evaporate, which resultr in the low 1.2-volt 
output. Typically, w.en the artery deprimes by, say, elevating the 
evaporator end , the output voltage first drops from 3.7 volts to an inter- 
mediate valve of 2.4 volts, and then again to 1.2 volts after the fillets 
evaporate. 

Priming tests were carried out by slowly lowering the evaporator 
end of the heat pipe in 0.1-inch steps until the output voltage from the 
hsat pipe indicated that a primed state was achieved. To check that the 
uninstrumented functional artery had also primed, 75 watts were applied, 
which results in a burnout if it had not. Of the four heat pipes, in two 
(X-1A and X-2C) the functional arteries had not yet primed. Further 
testing showed that a lower evaporator elevation by approximately 0.2 
inches was required to prime the functional arteries in these cases. 

These discrepencies can be attributed to the differences incurred during 
fabrication of the arteries, and the effect mentioned earlier (Section 5.3) 
of the reluctance of the liquid in the screen-walled portion of the artery 
to enter the foil-walled portion. In the case of heat pipe X-2C, the 
effect of thermistor cooling by liquid fillets ^"ound the thermistor 
appears to be greater than the other heat pipes to the extent that It 
is not clear whether the thermistor is detecting a primed state or fillets. 
Based on the above considerations, we have reconmended that heat pipes 
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X-1B and X-2D be designated as the primary ones for flight. The Pack-up 
heat pipes X-1A and X-2C have been recommended for the qualification 
vibration tests. 

5.6 Transient Tests with the Flight Power Profile 

The power profile to be used in the experiment is shown in 
Figure 5-4 beicw: 


100 WATTS 


J 1 1 I 1 1 

1 2 3 4 5 6 

TIME (MINUTES) 7V " roc - 

FIGURE 5-4. Power profile for the flight experiment. 

The duration of the zero-gravity period is expected to be 6 minutes. 

During the first half minute, 100 watts is applied to move excess liquid 
to the condenser end of the heat pipe. This amount of heat is sufficient 
to transfer 12% of the total amount of methanol. The zero-power period 
from t=.5 minutes to t=2 minutes is for priming. During the remainder 
of the zero-gravity period, 100 watts are again applied. If the functional 
artery has not sucessfully primed, a burnout will occur. 

Each heat pipe was tested repeatedly in the laboratory with 
th^ proposed power profile. To eliminate the contribution of puddle flow, 
which is not present in zero-gravity, the evaporator end of the heat pipe 
is elevated .5 inches for the final 100-watt period. During the transient 
tests, both temperature and thermistor voltage were recorded. The result 
of the transient test for heat pipe X-2D is shown in Figure 5-5. 

Initially the evaporator end is raised approximately 10 inches 
to ensure the arteries are in an unprimed state. At t»0 minutes, the 
heat pipe Is leveled and 100 watts is applied. The vapor temperature 
responds rapidly to the heat input and the thermistor voltage remains low. 
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indicating that the arteries are remaining unprimed. At t=.5 minutes 
the heat load is dropped to zero. The thermistor output voltage begins 
to climb as fillets form around the thermistor bead. Then abruptly at 
t *.8 minutes, the Instrumented artery primes. Film boiling occurs at 
the thermistor bead, as evidenced by the oscillations until t»1.3 minutes. 
At t»2 minutes, 100 watts is applied and the evaporator is elevated to 
.5 inches. At t=2.3 minutes the instrumented artery deprimes, as evidenc- 
ed by the step drop in thermistor voltage. The functional artery continues 
to operate, and the temperature difference between the vapor and condenser 
attains a steady-state value by t-4 minutes. At t= 4.8 minutes, fillets 
around the thermistor bead break away, and the thermistor voltage steps 
down again. The test terminates at t a 6 minutes. 

Several tests were run when the functional artery was 
Intentionally deprimed after t 2 minutes, and subsequently the evaporator 
temperature exceeds 200 ° F in less than one minute. 

5.7 Real-Time Power Profile Override 

We will have the opportunity to monitor data from the heat pipes 
during the flight and, if necessary, override the preprogrammed power 
profile. The data to be monitored is: 

o Evaporator temperature T^ measured by 
thermistor No. 1 (located 3 inches from 
the evaporator end). 

o Adiabatic temperature T 3 measured Dy thermistor 
No. 3 (located 16 inches from the evaporator 
end ). 

0 The voltage Vj H across the thermistor. 

The real-time control we require Is to override the programmed power 
profile with a zero-power command for any desired time Interval. The 
override will be exercised If certain conditions appear as the data is 
monitored. A list of the most probable anomalies Is given below: 

1. Vj H 5 0 (thermistor shorted). 

2 . V TH » 28 volts (thermistor lead broken). 
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3. 


4. 


5. 


6 . 


7. 


8 . 


9. 


10 . 


It takes longer than 2 minutes for Vj H to jump 
end indicate a primed state. 

A prime is indicated (Vy H >3 volts) at t = 0. 

An artery dump (drop in V TH ) is not indicated. 

Vy H is at some intermediate value or fluctuates in 
such a way that it is difficult to interpret. 

T-| - T 3 is not equal to zero at t = 0. 

T-| - T 3 does not drop to zero by t = 2 minutes. 

A burnout occurs after t * 2 minutes, (T-j -T 3 >20°F) 

T-j - T 3 may continuously increase after t * 2 minutes 
but not exceed 20° F. 


Of the above anomalies, we want to override the programmed power profile 
in three cases: Nos. 3, 8 and 9 above. The response to these anomalies 
are as follows: 

Anomaly Response 

3 If Vjh does not rise above 3.0 volts 

by t * 2 minutes, extend zero-power 
period to t = 3 minutes or unti 1 V TH 
does exceed 3.0 volts. 

8 If AT is greater than 3°F at t * 

2 minutes, extend zero-power period 
to t = 3 minutes or until AT is less 
than 3°F. 

9 If, after t * 2 minutes and before 

t * 4 minutes, AT exceeds 20° F, 
then command zero power until t * 

5 minutes. 

The first two responses above extend the priming period an extra minute 
If there is evidence that conditions for priming have not been achieved. 
The last response Is an attempt at a second try to prime the arteries if 
a burnout occurs. 
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7.0 NOMENCLATURE 


Section 2: 


A - 
B - 
k - 


n - 



< >.- 
' >c- 


< V 
( ),- 

( ) 2 - 


Area 

Multiplicative constant 

Boltzmann's constant 

Number of moles 

Partial pressure of gas 

Partial pressure of vapor 

Activation energy 

Gas constant 

Temperature 

Time 

Volume 

Refers to adiabatic section 

Refers to transition from parabolic to linear 
time dependence 

Refers to i th interval 
Refers to passivation period 
Refers to period beyond passivation 


Section 4: 

A - Cross-sectional vapor area in condenser 
c 

B - Constant 
c - Specific heat 
D - Diffusion coefficient 
g - Gravitational constant 
h^- Latent heat 
L c - Condenser length 
m - Mass 

M - Molecular weight 
N - Number of moles 

P - Pressure (also perimeter of heat-transfer area) 
Q - Heat 

R - Thermal resistance of condenser 
- Universal gas constant 
T - Temperature 
t - Time 
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U - Heat-transfer coefficient 

V - Total gas -blocked volume 

Z - Length of active portion of condenser 

c - Complex coefficient of t in exponent of assumed 

solution (Eq. (4-5)]. 

a,ui - Real and imaginary parts of c 

T p* T d " Constants (dimensions of time) 

t - Time lag 

p - Density 

( ) - Refers to adiabatic section 

a 

( ) c - Refers to condenser section 
( ) g - Refers to evaporator section 
(•) - Time rate 
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8.0 APPENDIX 

The Tables of Content of previous research reports follow. 
Research Reports Nos. 1 & 2 can also be found combi ned in Ref. [1-3]. 
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